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Abstract 

We report results from nearly simultaneous pentachromatic VLBI observa- 
tions towards a nearby GPS galaxy NGC 1052. The observations at 1.6 and 4.8 
GHz with VSOP, and at 2.3, 8.4, and 15.4 GHz with VLBA, provide linear reso- 
lutions of 0.1 pc. Convex spectra of a double-sided jet imply that synchrotron 
emission is obscured through foreground cold dense plasma, in terms of free-free 
absorption (FFA). We found a central condensation of the plasma which covers 
about 0.1 pc and 1 pc of the approaching and receding jets, respectively. A 
simple model with a geometrically thick plasma torus perpendicular to the jets 
is established to explain the asymmetric distribution of FFA opacities. 

Keywords: galaxies: active — galaxies: individual (NGC 1052, 0238 — 
084) — galaxies: nuclei — galaxies: jets — radio continuum: galaxies — 
techniques: high angular resolution 



1 Introduction 

GHz-Peaked Spectrum (GPS) sources show a convex radio spectrum peaked at GHz 
frequency, as are named (e.g. O'Dea 1998). The spectral shape indicates that a power- 
law spectrum of optically thin synchrotron radiation is affected by low-frequency cutoff. 
It is a controversial issue what is origin of the low-frequency cutoff; synchrotron self- 
absorption (SSA) or free-free absorption (FFA). 

Tight correlation between the peak frequency and the overall linear size (O'Dea & 
Baum 1997), and correlation between the equipartition size and the overall linear size 
(Snellen et al. 2000) indicates that SSA controls the the peak frequency, at which the 
equipartition size is derived by Scott & Readhead (1977). 
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Bicknell et al. (1997), on the other hand, pointed out that FFA through ionized gas 
surrounding radio lobes can also produce the correlation between the peak frequency 
and the overall size. Evidence for FFA was found via multi-frequency VLBI observa- 
tions in particular GPS sources; 1946+708 (Peck et al. 1999), OQ 208 (Kameno et al. 
2000; Xiang et al. 2002), 0108+388 (Marr et al. 2001), and NGC 1052 (Kameno et al. 
2001). 

We propose a clear method to discriminate the absorption mechanisms. In the 
case of SSA, the absorber is the synchrotron emitter itself such as jets or lobes. Sub- 
relativistic speeds of the emitter will cause significant Doppler shift of the peak fre- 
quency; approaching and receding jets will shows blue- and red-shifted spectra, respec- 
tively. Hence, we expect a higher peak frequency for an approaching jet. In the case 
of FFA, ambient cold dense plasma is an absorber whose speed is non-relativistic. We 
expect greater opacity towards receding jets, no matter how the absorber distributes, 
to produce a higher peak frequency. This picture was drawn towards some radio galax- 
ies such as 3C 84 (Vcrmeulen et al. 1994; Walker et al. 1994; Walker ct al. 2000), 
Cen A (Jones et al. 1996; Tingay & Murphy 2001), NGC 4261 (Jones et al. 2000; 
Jones et al. 2001), and NGC 6251 (Sudou et al. 2000). The opposite spectral be- 
havior between SSA and FFA allows us to discriminate the absorption mechanisms by 
measurements of the peak frequencies towards approaching and receding jets. Multi- 
frequency observations towards GPS sources with double-sided jets are relevant to this 
issue. 

We chose the nearest GPS source, NGC 1052, for this study. The redshift of this 
galaxy, z — 0.0049 (Knapp et al. 1978), corresponds to the distance of 20 Mpc, if we 
assume = 75 km s~^ and go = 0.5, and 1 milliarcsec (mas) corresponds to 0.1 pc. 
This object has parsec-scale twin jets in P. A. ~ 65° (Claussen et al. 1998), so that 
it is suitable for the discrimination of the absorbing mechanism. Their observations 
also revealed that H2O maser spots distribute along the western jet with a velocity 
gradient. 

Our previous trichromatic VLBA observations for the first time clarified the jet 
geometry and presence of a dense plasma torus that is responsible for FFA (Kameno 
et al. 2001). The eastern and western sides of the jet are approaching and receding 
with the velocity of p = 0.25 at the viewing angle of 50°. The western receding jet is 
obscured at lower frequencies, suggesting that a dense plasma torus covers the receding 
jet to produce FFA. Although the results were successful, three frequencies were not 
sufficient to obtain confident fit for the spectrum. 

In this paper, we report new pcntachromatic VLBI observations of NGC 1052, using 
the VSOP (Hirabayashi et al. 1998) and the VLBA. 

2 Observations and Results 

All pcntachromatic observations have been carried out within 18 days during July 
31 through August 31, 2001. We used the space- VLBI VSOP array at 1.6 and 4.8 
GHz, and VLBA at 2.3, 8.4, and 15.4 GHz. This combination provides comparable 
resolutions of ~ 1 mas through an order of frequency magnitude except at 2.3 GHz. 
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Figure 1: Pcntachromatic images of NGC 1052. Image performance is listed in table 1. 
Contours start at ±3 times the r.m.s. level and increase by factors of 2. Synthesized 
beams are shown at left bottom corner in each image. 



Observations at 2.3, 8.4, and 15.4 GHz using VLBA, coded BK084, were done 
on August 7, 2001. We used the dual-frequency 2.3/8.4 and the 15.4 GHz receiving 
systems, switching by 22-min scan cycles. Total integration time consists of 10 and 11 
scans for 2.3/8.4 and 15.4 GHz, respectively. Four 8-MHz baseband converters (BBCs) 
are allocated for 15.4 GHz, while 2 BBCs are used for 2.3 and 8.4 GHz, respectively. 
Fringes were detected with sufficient signal-to-noise ratios (SNRs) for all valid baselines 
and scans. 

Space- VLB! VSOP observations at 1.6 GHz (W513A) and 4.8 GHz (W513B) were 
executed on August 7 and July 31, 2001, respectively. The spacecraft HALCA was 
linked to two tracking stations during each observation. Green Bank and Tidbinbilla 
for W513A, and Green Bank and Usuda for W513B. Though fringes within ground 
array were strong enough, those in space baselines were marginal at 4.8 GHz. While 
visibilities with the tracking pass using Green Bank were available, wc missed them in 
the Usuda tracking pass. Fringes at 1.6 GHz were entirely detected within ground and 
space baselines. 

Fringe fitting, passband correction, and a priori amplitude calibrations were taken 
using the NRAO AIPS. Image synthesis and self calibration were processed with 
Difmap. Image performance is listed in table 1. 

Figure 1 shows images at five frequencies. Each image includes uncertainty in 
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J0241-08 at 4. 81 5 GHz 2001 Jul 31 




Right Ascension [mas] 




Frequency [GHz] Frequency [GHz] Frequency [GHz] 



Figure 2: Top: Total intensity map at 4.8 GHz convolved with an elliptical Gaussian 
whose FWHM is 1.9 x 1.4 mas in p.a.= 0°. Contours start at ±3cr level and increase 
by factors of 2, where a = 1.13 mJy beam~^. Approaching and receding components 
are marked as and R^, where k = 1, 2, and 3. Bottom: Spectra of components. 



absolute position through self-calibration process. To register the images at different 
frequencies, we picked up 6 distinct components; Al, A2, and A3 in the approaching 
jet, and Rl, R2, and R3 in the receding jet. We then measured positions {^^,i]k) of 
them with respect to a tentative origin of each image. Here k and u stand for indices 
of components and frequency, respectively. We then derived relative offsets {6^,'^, Sr]") 
to minimize the positional residuals xl defined as 

where a^^ and o"^^ are standard positional errors of component k at frequency u, and 
1^0 is the frequency of the reference image. We take uq = 4.8 GHz, where the frequency 
is middle of our observing range, and where the synthesized beam is the narrowest. 
Standard errors of the registration are listed in table 2. 



k 
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Table 1: Image performance 



Frequency 




Beam 




Image r.m.s. 


Telescopes 


(GHz) 


Major (mas) 


Minor (mas) 


P.A. (deg.) 


(mJy beam~^) 




1.6 


1.75 


1.39 


31.3 


1.183 


HALCA, ATCA, and VLBA 


2.3 


6.34 


2.58 


-5.8 


0.511 


VLBA 


4.8 


0.78 


0.67 


-78.7 


1.473 


HALCA and VLBA 


8.4 


1.86 


0.77 


-2.2 


0.444 


VLBA 


15.4 


1.09 


0.42 


-7.8 


0.432 


VLBA 



Table 2: Accuracy of image registration 



Frequency 


R. A. 


Dec. 


(GHz) 


(mas) 


(mas) 


15.4 


0.30 


0.29 


8.4 


0.21 


0.31 


4.8 






2.3 


2.80 


0.75 


1.6 


0.94 


0.43 



3 Discussion 

First of all, we aim at clarifying the absorbing mechanism. Figure 2 shows continuum 
spectra of 6 components. Receding components Rl, R2, and R3 show sharp cutoff 
at low-frequency ends. The spectral index of R2 between 1.6 and 2.3 GHz is afl — 
4.76 ± 0.48. Here we define the spectral index a as oc z/"*"", where is the fiux 
density and f is the frequency. The spectral index exceeds maximum attainable spectral 
index of 2.5 by SSA. For components R2 and R3 j > 2.56 and af^ = 2.43 ± 0.38, 
respectively, are close to the SSA hmit. We evaluated the spectral indices between 1.6 
and 4.8 GHz for these components because flux densities at 2.3 GHz is too weak to 
detect. The brightest approaching component, A3, shows af^ = 3.28 ± 0.05. These 
steeply rising spectra indicate signiflcant presence of FFA, which can produce a > 2.5. 

Component Al shows a single power-law spectrum with a — —1.1, whose spectral 
peak should be at < 1.6 GHz, while the spectral peak of its counter component Rl 
is between 2.3 and 4.8 GHz. Components A2 and A3 mark the spectral peaks at 
~ 2.3 and ~ 8.4 GHz, respectively, whose counter components R2 and R3 show the 
spectral peaks at ~ 8.4 and > 15.4 GHz. Every approaching component shows a lower 
peak frequency than each counterpart does. As is mentioned in the introduction, this 
manner clearly indicates that the receding jet is obscured by external absorber, say, 
FFA. 

To illustrate spatial distribution of opacity, we applied spectral fit for the FFA 
model 

5, = 5oi^"Oexp(-roz/-'-i), (2) 
pixel by pixel. Here, Si, is the observed fiux density at the frequency u in GHz, ao is 
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NGC1052 




Figure 3: Top: Spatial distribution of FFA opacity tq in grey scale overlaid to the total 
intensity map at 4.8 GHz in contour. Bottom: Profile of the FFA opacity along the 
jet. Errors are derived from the nonlinear spectral fit using equation 2. 

the intrinsic spectral index, 5*0 is the absorption-corrected flux density at 1 GHz, and 
To is the FFA coefficient. We fixed ckq to —1.1, which is obtained from the spectral 
index of component Al. Figure 3 shows spatial distribution of tq and its profile along 
the jet. The maximum value of tq ~ 300 is consistent with Kameno et al. (2001). The 
profile tails over ~ 1 pc towards the western receding jet, while it rapidly decreases 
towards the eastern approaching jet. This asymmetric profile is again consistent with 
Kameno et al. (2001), with more detailed structure at the hands of high resolution by 
VSOP. 

Since the opacity coefficient, tq, is related to the electron density, Ue, and temper- 
ature, Tg, along line of sight as 

To = 0.46 / nlTe^L, (3) 

JLOS 

it will be an indicator of cold dense plasma. Here, the integration stands for path 
length along the line of sight in unit of pc. When the aborber is homogeneous, tq is 
proportional to the path length through the ionized region. A simple model to explain 
the opacity profile is shown in figure 4 with a geometrically thick torus perpendicular 
to the jets. The inner surface of the torus is photo-ionized by illumination of ionizing 
photons from the nucleus, while outer region remains in neutral gas or dusts. The path 
length through the ionized region can produce the asymmetric profile of the opacity, 
when we look from the viewing angle of 50°. 

We can estimate physical condition of the absorber, if we assume the torus model. 
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Figure 4: Schematic diagram of NGC 1052. Double-sided jet is inclined by 50° with 
respect to the line of sight. A geometrically thick molecular torus is perpendicular to 
the jet, whose inner surface is illuminated by ionizing photons from the nucleus. The 
ionized region obscures radio emission from the jet, in terms of free-free absorption. 
Long path length through the ionized region towards the nucleus results in a large 
opacity, as is shown in figure 3 

The maximum value of tq ~ 300 towards nucleus gives 



when the path length is the same the width of the absorbed area along the receding 
jet, say, L ^ 1 pc. The ionizing condition of > 10^ K gives the lower limit of 
He > 2.5 X 10^ cm~^. The upper limit of Tg can be given by absence of significant 
free-free emission below the brightness temperature of Tb < 4.9 x 10^ K at 15.4 GHz 
or Tb < 5.4 X 10^ K at 8.4 GHz. Since the torus is optically thick, in terms of 
bremsstrahlung, we have Tg ~ Tb < 4.9 x 10^ K and then Tig < 2.6 x 10^ cm~^. This 
evaluation indicates that the FFA absorber consists of cold (10^ < Tg < 4.9 x 10^ K) 
and dense (2.5 x 10^ < rig < 2.6 x 10^ cm plasma. The average value of electron 
density derives the electron column density across the plasma torus (L ~ 1 pc). 



This electron column density is comparable to the atomic column density of ~ 10^'^ 
cm~^ measured by ROSAT and ASCA X-ray observations (Guainazzi & Antonelli 
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(4) 
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(5) 



1999). 
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We mention that the location of water maser spots (Claussen et al. 1998) coincides 
with where the FFA opacity is large. The torus model explains results of the water 
masers; the H2O molecules survive inside the torus, shielded against exposure by the 
nucleus. The seed photon from the receding jet are amplified by excited molecules in 
the torus. Amplification does not work towards the approaching jet, because the torus 
is not in front of the jet, the seed photon supplier. Absence of masers towards the 
eastern jet can be understood by this model. A detailed monitoring program would 
verify this model by detection of orbital motion of maser spots within the torus. 
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